Searching for a Cosmological Preferred Axis: 
Union2 Data Analysis and Comparison with Other Probes. 
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We review, compare and extend recent studies searching for evidence for a preferred cosmological 
axis. We start from the Union2 Snia dataset and use the hemisphere comparison method to search 
for a preferred axis in the data. We find that the hemisphere of maximum accelerating expansion 
rate is in the direction (/, b) = (309°_3o , 18°_]^qo) (r^om = 0.19) while the hemisphere of minimum 

acceleration is in the opposite direction {l,b) — (129°l3o ,— 18°lJ5o) (^Om = 0.30). The level of 
anisotropy is described by the normalized difference of the best fit values of Qom between the two 
hemispheres in the context of ACDM fits. We find a maximum anisotropy level in the Union2 data 



of 



0.43 ±0.06. Such a level does not necessarily correspond to statistically significant 



anisotropy because it is reproduced by about 30% of simulated isotropic data mimicking the best fit 
Union2 dataset. However, when combined with the axes directions of other cosmological observations 
(bulk velocity flow axis, three axes of CMB low multipole moments and quasar optical polarization 
alignment axis), the statistical evidence for a cosmological anisotropy increases dramatically. We 
estimate the probability that the above independent six axes directions would be so close in the sky 
to be less than 1%. Thus either the relative coincidence of these six axes is a very large statistical 
fluctuation or there is an underlying physical or systematic reason that leads to their correlation. 



PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



1. INTRODUCTION 

A large and diverse variety of cosmological observa- 
tions during the past twenty years have established a 
standard cosmological model (ACDM ) which is based on 
the cosmological principle (homogeneity, isotropy, gen- 
eral relativity, cold dark matter + baryonic matter), flat- 
ness of space, the existence of a cosmological constant 
and gaussian scale invariant matter perturbations gen- 
erated during inflation. This cosmological model makes 
clear and well defined predictions which have withstood 
the continuous and rapid improvement of cosmological 
observational tests. Prominent successes of the standard 
cosmological model include the following: 

• The Cosmic Microwave Background (CMB) angu- 
lar power spectrum of perturbations [1 is overall 
in excellent agreement with the predictions of the 
standard model. However, a few issues related to 
the orientation and magnitude of low multipole mo- 
ments (CMB anomalies) constitute remaining puz- 
zles for the standard model |2j-t9j. 

• The statistics of the CMB temperature perturba- 
tion maps [10 are consistent with the prediction of 
gaussianity of the standard model. 

• Observations of the recent accelerating expansion 
history of the universe [11 are consistent with the 
existence of a cosmological constant. Despite of the 
continuously improved data no need has appeared 
for more complicated models based on dynamical 
dark energy or modified gravity. The likelihood 
of the cosmological constant vs more complicated 
models has been continuously increasing during the 



past decade [l2]. 

• Observations of large scale structure are in good 
agreement with ACDM [13 (basic statistics of 
galaxies [15 , halo power spectrum [16 ). 

Despite of the above major successes the standard model 
is challenged by a few puzzling large scale cosmological 
observations [M] which may hint towards required mod- 
ifications of the model. These challenges of ACDM may 
be summarized as follows: 

1. Large Scale Velocity Flows: ACDM predicts 
significantly smaller amplitude and scale of flows 
than what observations indicate. It has been 
found that the dipole moment (bulk flow) of a 
combined peculiar velocity sample extends p/7j on 
scales up to 100h~^Mpc {z < 0.03) with ampli- 
tude larger than 400km/sec. The direction of 
the flow has been found consistently to be ap- 
proximately in the direction / c^ 282°, b c^ 6° . 
Other independent studies have also found large 
bulk velocity flows on similar directions |18j on 
scales of about 100h~^Mpc or larger [19]. The ex- 
pected rms bulk flow in the context of ACDM nor- 
malized with WMAP5 (l^om,cr8) = (0.258,0.796) 
on scales larger than 50h~^Mpc is approximately 
110km/ sec. The probability that a flow of magni- 
tude larger than AOOkm/sec is realized in the con- 
text of the above ACDM normalization (on scales 
larger than bOh~^Mpc) is less than 1%. A possi- 
ble connection of such large scale velocity flows and 
cosmic acceleration may be found in Ref. [ 20] . 

2. Alignment of low multipoles in the CMB 
angular power spectrum: The normals to the 



octopole and quadrupole planes are aligned with 
the direction of the cosmological dipole at a level 
inconsistent with Gaussian random, statistically 
isotropic skies at 99.7% [3 . The corresponding 
directions are: octopole plane normal (/, 6) = 
(308°, 63°) [a in], quadrupole plane normal (/, 6) = 
(240°, 63°) 0122], CMB dipole moment {l,b) = 
(264°, 48°) [23]. A related effect has also been re- 
cently observed [24] by considering the tempera- 
ture profile of 'rings' in the WMAP temperature 
fluctuation maps. It was found that there is a ring 
with anomalously low mean temperature fluctua- 
tion with axis in the direction (/,6) = (276°,—!°) 
which is relatively close to the above directions 
(particularly that corresponding to the bulk veloc- 
ity flows). 

3. Large scale alignment in the QSO optical po- 
larization data: Quasar polarization vectors are 
not randomly oriented over the sky with a probabil- 
ity often in excess of 99.9%. The alignment effect 
seems to be prominent along a particular axis in 
the direction {l,b) = (267°, 69°) [25^. 

4. Profiles of Cluster Haloes: ACDM predicts 
shallow low concentration and density profiles in 
contrast to observations which indicate denser high 
concentration cluster haloes [26l [271 • 



5. Missing power on the low / multipoles of the 
CMB angular power spectrum which leads to 
a vanishing correlation function C{0) on angular 
scales larger than 60° 



In addition to the above large scale effects there are is- 
sues on galactic scales (missing satellites problem [29-31 
and the cusp/core nature of the central density profiles 
of dwarf galaxies [32 ). In what follows we focus only on 
the kind of anomalies that appear to be related with the 
possible existence of a preferred axis. The rest are only 
mentioned for completeness but we can not exclude the 
possibility that a single mechanism could be responsible 
for seemingly unrelated anomalies. After all, before in- 
flation, it would be hard to imagine that the same mech- 
anism could be responsible for both the generation of 
primordial fluctuations on all scales and for the horizon 
problem (two seemingly unrelated issues). 

Three of the above five large scale puzzles are large 
scale effects related to preferred cosmological directions 
(CMB multipole alignments, QSO polarization alignment 
and large scale bulk flows) which appear to be not far 
from each other. Their direction is approximately nor- 
mal to the axis of the ecliptic poles (/, 6) = (96°, 30°) and 
lies close to the ecliptic plane and the equinoxes. This 
coincidence has triggered investigations for possible sys- 
tematic effects related to the CMB preferred axis but no 
significant such effects have been found [3 . 

Thus, unless there is a hidden common systematic [33], 
the existence of a cosmological preferred axis may be at- 



tributed to physical effects. An incomplete list of these 
effects is the following: 

• An anisotropic dark energy equation of state [34]- 
[36] due perhaps to the existence of vector fields 

[371 Eg. 

• Dark Energy and/or Dark matter perturbations on 
scales comparable to the horizon scale [39^ .40^ . For 
example an off center observer in a IGpc void would 
experience the existence of a preferred cosmologi- 
cal axis through the Lematre-Tolman-Bondi metric 
[41H45]. 

• Deviations from the isotropic cosmic expansion rate 
induced by a fundamental violation of the cosmo- 
logical principle eg through a multiply connected 
non-trivial cosmic topology [46l |47] , rotating uni- 
verse coupled to an anisotropic scalar field [48] . 
non-commutative geometry [49] or simply a fun- 
damental anisotropic curvature f50^. 

• Statistically anisotropic primordial 
perturbations [5TH54] . For example, inflation- 
ary perturbations induced by vector fields [55H58] . 
Note however that inflationary models with vector 
fields usually suffer from instabilities due to the 
existence of ghosts. [59] 

• The existence of a large scale primordial magnetic 
field [60- 62 . Evidence for such a magnetic field has 
recently been found in CMB maps [63 . 

In view of the convergence of the directions of the pre- 
ferred axes mentioned above, it is important to inves- 
tigate other cosmological datasets for hints of existence 
of preferred axes. Type la supernovae used as standard 
candles to investigate the recent expansion history of the 
universe constitute one such class of datasets. 

Most previous studies searching for anisotropics in 
Snia datasets have found no statistically significant evi- 
dence for anisotropics [64 £9, . The latest such study [65] 
used the latest and largest dataset available (Union2 [11] 
consisting of 557 SnIa), to derive the angular covariance 
function of the standard candle magnitude fluctuations 
searching for angular scales where the covariance func- 
tion deviates from in a statistically significant manner. 
No such angular scale was found. This is a useful and am- 
bitious approach which aims at identifying not only the 
existence of a possible anisotropy but also its detailed an- 
gular scale dependence despite of relatively small number 
of data in most angular scales considered. 

An alternative approach is found in [70 where a sta- 
tistically significant preferred axis was found using the 
hemisphere comparison method. This method, adopted 
also in the present study (see next section), amounts to 
fitting the ACDM parametrization for (^om on several 
pairs of opposite hemispheres and comparing the maxi- 
mally asymmetric pair with the corresponding maximally 
asymmetric pair of isotropized similar datasets. The ad- 
vantage of this method is that it optimizes the statistics 



since there is a large number of Snia in each hemisphere. 
This is achieved at the cost of loosing all information 
about the detailed structure of the anisotropy. 

The method was applied to four SnIa datasets [70^. 
The most prominent axis of maximal hemispheric asym- 
metry (/,6) = (123°, 27°) was found to be close to the 
axis of the equatorial poles (/, 6) = (96°, 30°). This align- 
ment between maximally asymmetric hemispheric Hub- 
ble diagrams and the equatorial frame was attributed to 
a systematic error by the authors of Ref. [70 . A second 
direction of maximum asymmetry (/,6) = (235°, 15°) in 
one of the datasets (the Gold04 [71]) considered in Ref. 
[70] was closer to the preferred axes of other observa- 
tions discussed above but it was also discarded because 
the maximum asymmetry between the hemispheres co- 
incided with the maximal asymmetry in the number of 
degrees of freedom. 

In this study we apply a variation of the hemisphere 
comparison method to the Union2 dataset. Our goal is to 
identify the direction of the axis of maximal asymmetry 
for the Union2 dataset [11] (directions to the SnIa pro- 
vided in Ref. [65 ) and evaluate its statistical significance 
by comparing with a large number of similar Monte Carlo 
isotropized datasets. We also compare the obtained di- 
rection of maximal asymmetry axis with the directions 
of preferred axes obtained with the other cosmological 
observations discussed above. In particular we find the 
probability that these axes would have the observed an- 
gular separation if they were uncorrelated. 

The structure of this paper is the following: In the next 
section we describe in more detail the hemisphere com- 
parison method and we apply it to the Union2 dataset. In 
section 3 we compare the direction of the axis of maximal 
asymmetry with the directions of preferred axes from the 
other observations discussed above. Finally in section 4 
we conclude summarize and discuss future prospects of 
this work. 



rial coordinates (right ascension and declination) of each 
SnIa. It is straightforward to convert these coordinates 
to galactic coordinates or to usual spherical coordinates 
(^, ^) in the equatorial or galactic systems [75^. 

This dataset may be analyzed in the usual manner by 
applying the maximum likelihood method. The apparent 
magnitude m{z) is related to the Hubble- free luminosity 
distance Dl{z) through 

mth{z, l^om) = M(M, Ho) + 5logio{DL{z)) (2.1) 

where in a fiat cosmological model 

DUz) = {l + z) Tdz' ^° (2.2) 

Jo H[z';\lom) 

is the Hubble-free luminosity distance assumed here to 
be parameterized by ACDM i.e. 



H{zf = ifo'["om(l + Zf + (1 - l^Om) 



(2.3) 



Also M is the magnitude zero point offset which depends 
on the absolute magnitude M and on the present Hubble 
parameter Hq as 



M = M + 5/o^io(^^)H 
Mpc 

= M -blogioh ^42.38 



25 



(2.4) 



The parameter M is the absolute magnitude which is 
assumed to be constant. 

The data points of the Union2 dataset are given in 
terms of the distance modulus 



f^obsizi) = rriobsizi) - M 



(2.5) 




2. Evaluate Best Fit Q^ in each Hemisphere 



2. MAXIMUM ASYMMETRY AXIS OF THE 
UNION2 DATASET 

The Union2 SnIa dataset [11 is a complilation con- 
sisting of 557 SNe la covering the redshift range z = 
[0.015, 1.4]. It extends the Union dataset [72j by adding 
new SnIa data at low and intermediate redshifts discov- 
ered by the CfA3 [73 and SDSS-II Supernova Search 
[74] . respectively. It also includes six new SnIa discov- 
ered by the Hubble Space Telescope at high z. Here we 
use the directions to the SnIa provided in Ref. [65 . The 
angular distribution of the Union 2 dataset in galactic co- 
ordinates, is shown in Figure 1. The color of each point 
provides information about the redshift. 

The Union2 data along with directions as presented 
in Ref. [£5j include the SnIa name, the redshift in the 
CMB rest frame, the distance modulus and its uncertain- 
ties (which include both the observational and the intrin- 
sic magnitude scatter). They also include the equato- 




Union2 Data 

Galactic Coordinates 

(view of sphere from opposite directions 



4. Repeat with several random axes 
and find AQ,,,^,^ 



Fig. 1: The Union2 data in galactic coordinates and the Hemi- 
sphere Comparison Method. The color of each datapoint pro- 
vides information about its redshift according to the legend 
on the right. Two opposite hemispheres are shown. The view- 
points are along the axis of maximum asymmetry discussed 
below. 



The theoretical model parameter (r^om) is determined by 
minimizing 



X^(^Om,Mo) 



Y^ {f^obsj^i) — Mt/i(^i7^0m7Mo)) /^ a.\ 



= 1 



a 



ji I 



where a^ ^ are the distance modulus uncertainties which 
include both the observational and the intrinsic magni- 
tude scatter. These uncertainties are assumed to be gaus- 
sian and uncorrelated (we assume a diagonal covariance 
matrix and ignore systematics). The theoretical distance 
modulus is defined as 

(2.7) 
where 



/io = 42.38 — blogioh 



(2.8) 



As a test we have checked that our full sky analysis 
reproduces the results of Ref. [11 for the cases of no 
systematics. For example we obtain a full sky best fit 
value l^om = 0.27 and we reproduced the contour Fig. 
10a of Ref. [IT]. 

The hemisphere comparison method implemented in 
our study involves the following steps (see also Fig. 1): 



1. Generate a random direction 



'f^rnd = (cos (j) \/l — U^ , siu (j) V 1 — U^ ^u) 



(2.9) 



where (j) G [0, 27r) and u ^ ["l^ 1] are random num- 
bers with uniform probability distribution. 



Anisotropies for Random Axes (UnionZ Data) ^ 
View from above Maximum Asymmetry Axis 
Galactic Coordinates 



Minimum Acceleration: ^^m.. 
(l,b)=(129'',-18°) ^ 





I An/Q= 

I 0.43 



IVIaximum Acceleration Direction: 



Anl!P 

-0.43 



Fig. 2: The directions of the random axes considered are 
shown as dots on the unit sphere colored according to the 

sign and magnitude of the anisotropy level ( ^ °"^ J . The 

hemisphere shown on the left (right) is the one correspond- 
ing to maximum (minimum) acceleration. The corresponding 
best fit values of Qom are r^om = 0.19 ((/,6) = (309°, 18°)) 
and Qom = 0.30 ((/, h) = (129°, -18°)). 



2. Split the dataset under consideration into two sub- 
sets according to the sign of the inner product 
^rnd ' ^dat whcrc r^at IS a uuit vector describing 
the direction of each Snia in the dataset. Thus one 
subset corresponds to the hemisphere in the direc- 
tion of the random vector (defined as 'up') while 
the other subset corresponds to the opposite hemi- 
sphere (defined as 'down'). 

3. Find the best fit values on l^om in each hemisphere 
(^Om,n and ^om,d)- Use these values to obtain the 
anisotropy level quantified through the normalized 
difference 



An 



Om 



n 



Om/ 



n 



Om,d 



n^ 



Om 



^Om,n + ^On 



(2.10) 



4. Repeat for 400 random directions Trnd and find the 
maximum standardized difference for the Union2 
data ( ^^0"^ rn,aa: \ ^ ^g g^^gQ obtalu the corre- 






sponding direction of maximum anisotropy. 

Instead of the above algorithm involving random direc- 
tions of the axes, we could have implemented a uniform 
coverage of the sphere, utilizing equal area pixels along 
the lines of Healpix [76]. In that case each axis would 
have a fixed direction in the center of each equal area 
pixel. However, we find the use of axes with random 
directions simpler to implement in practice without any 
specific disadvantage compared to the Healpix approach. 
In order to maximize efficiency, the number of axes 
should be approximately equal to the number of data 
points (SnIa) on each hemisphere. The reason for this is 
that changing the direction of an axis, does not change 



the corresponding 



A^o 



^Or, 



until a data point is crossed 



by the corresponding equator line. Such a crossing is ex- 
pected to occur when the direction of an axis changes by 



Maximum Acceleration Directions 
for Different Redshift Cutoffs 



IN/laximum Acceleration Direction: 
(l,b)=(341'',-17'>) (Zmax=0.2) 





IVIaximum Acceleration Direction: 
(l,b)=(309°,18'') (full sample) 



Fig. 3: The directions of maximum acceleration corresponding 
to the redshift cutoffs of Table 1 . 



approximately the mean angular separation between data 
points. Thus, using more axes than the number of data 
points in a hemisphere does not improve the accuracy of 
the determination of the maximum anisotropy direction. 
Given that the number of datapoints per hemisphere for 
the Union2 dataset is about 280 it becomes clear that us- 
ing 280 test axes is close to the optimal number of axes 
to use. We have used used 400 axes in our analysis, well 
above the value of 280. 

In order to derive the la error in the maximum 
anisotropy direction, we first obtain the la error ctaq 



associated with 



O'AQ 



A^o 



- . This is of the form 






a 



Om.n 



^ 



= 0.06 



Om,d 



(2.11) 



Notice that this is the error due to the uncertainties of 
the supernova distance moduli propagated to the best fit 
ftorn on each hemisphere and thus to — ^'"^"^ . We then 
identify all the test axes that correspond to an anisotropy 
level within la from the maximum anisotropy level ie 



AOn 



A^o 



^ — ^ - ± ctaq. These axes apparently cover 
an angular region corresponding to the la range of the 
maximum anisotropy direction. Using a run with 400 test 
axes we find the la angular region for the hemisphere of 
maximum acceleration is in the direction / = 309° ^30 , 

b = 18°l]^Jo (best fit l^om = 0.19) while the hemisphere 
with the minimum acceleration is in the opposite direc- 
tion (/,6) = (129°1^3o,-18°llfo) (best fit Qom = 0.30). 
The corresponding maximum anisotropy level is 



AU 



Om max 



a 



Om 



U2 



0.43 ± 0.06 



(2.12) 



Our results are shown in Fig. 2 where we present the 
directions of the random axes considered, as dots on the 
unit sphere colored according to the sign and magnitude 



U2 



of the anisotropy level ( ^ °"^ J . The view point is on 

top of the hemisphere of maximum (left) and minimum 
(right) acceleration. 

In an effort to identify possible redshift dependence 
of the above anisotropy, we have implemented a redshift 
tomography of the Union2 data and have identified the 
maximum anisotropy directions (with their errors) for the 
following redshift ranges: 0-0.2, 0-0.4, 0-0.6, 0-0.8, 0-1.0, 
0-1.2. 0-1.4 (see Table 1 and Fig. 3). We have found 
that with the exception of the 0-0.2 redshift range whose 
maximum anisotropy direction is about 40° away from 
the maximum anisotropy direction of the full dataset, all 
the other redshift ranges have an anisotropy direction 
which is within about 20° from the anisotropy direction 
of the full Union2 dataset. 

We next wish to addres s the question whether the 
maximum anisotropy level (2.12) for the Union2 data is 



consistent with statistical isotropy. In order to address 
this question we have constructed simulated isotropic 



datasets by replacing the v'^ distance modulus of the 
Union2 dataset by a random number with a gaussian 
distribution with mean and standard deviation deter- 



mined by the best fit value of iith{zi^ ^om, Mo) (eq. (2.7) 
with the best fit values l^om = 0.27, /io = 43.16) and 
by cr^ i of the corresponding Union2 data-point respec- 
tively. We then compare a simulated isotropic dataset 
with the real Union2 dataset by spliting each dataset 
into hemisphere pairs using 10 random directions. This 
part of our analysis is not aimed at identifying the max- 
imum anisotropy direction neither at comparing with 
the result of the search in N=400 directions. Instead it 
only aims at comparing the real data with the isotropic 
simulated data with respect to the level of anisotropy. 
In the context of this comparison it is not important 
to identify the level of absolute maximum anisotropy. 
What is more important is to repeat the comparison 
a relatively large number of times (40 in our case) in 
order to have acceptable statistics. Given the limita- 
tions of computing time we had to reduce the number of 
axes directions (10) in order to increase the number of 
Monte Carlo experiments performed. Clearly the level 
of anisotropy identified in this case ( ^ "^ ^ 0.29) is sig- 
nificantly smaller compared to the case of 400 axes di- 
rections (^p^ — 0.43) but this is not important for our 
purposes which do not include in this case the identi- 
fication of the maximum anisotropy. Next we find the 

/ xC/2 

maximum levels of anisotropy ( ^"^ ^^^ ) (Union2) 

and ( — ^ ^^^ J (Isotropic) and compare them. We re- 
peat this comparison experiment 40 times with different 
simulated isotropic data and axes each time. We found 
the following results: 

• In about 1/3 of the numerical experiments (14 

times out of 40) ( ^^Om ^a:, \ y f A^Om max ] ^g_ 

the anisotropy level was larger in the isotropic sim- 
ulated data. In the rest 2/3 of the numerical exper- 
iments the anisotropy level was larger in the Union2 
data. This is a clear indication that the anisotropy 
level found in the Union2 data is consistent with 
statistical isotropy. 

• The mean and standard deviation of the maximum 
anisotropy levels in each case are 



A^ 



Om max 



a 



Om 



U2 



for the Union2 data and 



a 



Om 



0.29 ±0.05 



0.24 ±0.07 



(2.13) 



(2.14) 



for the simulated isotropic data. The error region 
described by equation (2.13) is not the error associ- 



ated with the uncertainties of the supernova mag- 
nitudes as are the error regions of the 4^^ and 5*^ 



Table 1: Directions of maximum anisotropy for several redshift ranges of the Union2 data (see also Fig. 3). The level of 
maximum anisotropy for a typical isotropic simulated dataset is also shown in the 5*^ column. The asymmetry of errors is 
largely due to the non-uniform distribution of the Snia on the sky. 



Redshift Range 


/ 


b 




/^^ \SIM 

( A^om max ] 
V ^Om J 


0-0.2 


341°1^22 


-17°lf 


2.08 it 0.22 




4.28 ±0.22 


0-0.4 


301°!^' 


10+21 

-J- -15 


1.81 ±0.27 




1.23 ±0.15 


0-0.6 


sorttl 


-4°ir6 


0.6 ±0.1 




0.48 ±0.08 


0-0.8 


327°1^2? 


srltg 


0.46 ±0.07 




0.36 ±0.06 


0-1.0 


301°lg^ 


_40+31 


0.45 ±0.07 




0.35 ±0.06 


0-1.2 


310°+4 


16°+1? 


0.43 ±0.07 




0.35 ±0.06 


0-1.4 (all data) 


309°+f 


-100+11 
-•-» -10 


0.43 ±0.06 




0.36 ±0.06 



columns of Table 1. It corresponds to the range 
of the anisotropy level obtained when using 10 test 
axes to find the maximum anisotropy. Given the 
relatively small number of axes directions (10) con- 
sidered in this part of the analysis, there is a signif- 
icant variation of the level of maximum anisotropy 
identified in each run . Th is variation is described 
by the result of eq. (2.13). In this case what we 



call 'the error' corresponds to the range of values 
— ^T'"^''"' that is expected to be obtained in about 
68% of the trials when our approach is implemented 
using 10 axes with random directions. In view of 
the fact that this variation is large enough to over- 



Distribution of 



Isotropic Data can Reproduce 

the Asymmetry Level of 

Union2 Data 



Maximum ^^m 

Q 
after 10 trial axes 




Fig. 4: The distribution of the anisotropy levels ( ^^^ ^^^ j 

in the Union2 (front histogram) and simulated Isotropic 
datasets (back) using 10 random directions in each dataset. 
They show significant overlap which is a sign of consistency 
of the Union2 data with statistical isotropy. 



lap significantly with the corresponding results ob- 
tained with the isotropic Monte Carlo data, it be- 
comes clear that the anisotropy of the real data is 
consistent with statistical isotropy. It is clear from 
eqs. (2.13) and (2.14) that there is a clear overlap at 



the la level which also implies that the maximum 
anisotropy level of the Union2 data is consistent 
with statistical isotropy. 

• The histograms indicating the distribution of 

( ^^0"^ ^^^ ) in each case are shown in Fig. 3 and 

they clearly show a significant overlap confirming 
also the consistency of the Union2 data with sta- 
tistical isotropy. 

Thus we have identified a direction of maximum 
anisotropy in the Union2 data and the level of this 
anisotropy is larger than about 70% of isotropic simu- 
lated datasets. However this level is clearly not enough 
to indicate inconsistency with statistical isotropy. 

In the next section we proceed to compare the direction 
of the identified maximum anisotropy with correspond- 
ing directions obtained with different, independent cos- 
mological observations. The potential consistency among 
these independent anisotropy directions can dramatically 
increase the statistical significance of each one of them. 



3. CORRELATIONS WITH PREFERRED AXES 
FROM OTHER OBSERVATIONS 

As discussed in the Introduction, there is a range of in- 
dependent cosmological observations which indicate the 
existence of anisotropy axes. The consistency of these di- 
rections may be interpreted as a hint of the existence of 
an underlying physical or systematic cause which is com- 
mon in all of these apparently independent axes. These 
cosmological observations along with their preferred di- 
rections and the corresponding references are summa- 
rized in Table 2. In what follows we focus on the preferred 



Table 2: Directions of Preferred axes from different cosmo- 
logical observations. The precise directions may vary by a few 
degrees across the hterature but our results are insensitive to 
such small variations. 



Cosmological Obs. 


/ 


b 


Reference 


Snia Union2 


309° 


18° 


This Study 


CMB Dipole 


264° 


48° 


m 


Velocity Flows 


282° 


6° 


mm 


Quasar Alignment 


267° 


69° 


m 


CMB Octopole 


308° 


63° 


m 


CMB Quadrupole 


240° 


63° 


[22] ED 


Mean 


278° ± 26° 


45° =b 27° 


- 



it). It is straightforward to estimate the probability that 
six random points would lie is such a small region on 
a hemisphere. To obtain this estimate we evaluate the 
mean value of the inner product between all pairs of unit 
vectors corresponding to the preferred directions of Table 
2 and compare with the corresponding mean value of six 
random directions on a hemisphere. Thus we evaluate 



< \cosOij\ >=< 



Vi -rA >= 



N 



^^ N(N- 



1) 



(3.1) 



where in our case N = 6 and we take the absolute value 
because we are ignoring the directionality of the axes. 
We thus apply eq. (3.1) to both the real data of Table 
2 and to 1000 realizations of six random points on the 
sphere obtained using eq. (2.9). For the real data we 
find 



axes and ignore the information about the directionality 
of each axis. Had we considered also information about 
the directionality of axes the likelihood of the observed 
coincidence would be even smaller given that the direc- 
tions of the bulk velocity flow, faster accelerating expan- 
sion and CMB dipole all appear to be towards the North 
Galactic Hemisphere. 

The six axes corresponding to the observations of Table 
2 are shown in Fig. 4 in galactic coordinates. Clearly 
all coordinates are located in a relatively small part of 
the North Galactic Hemisphere (less than a quarter of 
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Fig. 5: The coordinates of the preferred axes of Table 2 are all 
located in a region less than a quarter of the North Galactic 
Hemisphere (left). The south galactic hemisphere (right) is 
also shown for completeness. The bulk flow direction is also 
visible in the south galactic hemisphere because it is close to 
the equator. The mean direction obtained in Table 2 with 
coordinates (/,6) = (278°, 44°) is also shown. 



< \cosOij\ >=0.72 (3.2) 

while from the Monte Carlo uncorrelated data we obtain 



< \cosOij\ >= 0.5 ±0.072 



(3.3) 



Clearly, the value of < \cos6ij \ > of the real data is about 
3cr away from the expected value if there were no corre- 
lation among the axes of Table 2. This is also seen in 
Fig. 5 which shows a histogram of the distribution of 
< \cosOij\ > as obtained from the Monte Carlo data su- 
perposed with the value obtained from the real data (eq. 
( |3.2| ). Clearly, less than 1% (0.8%) of the simulated data 
exceed the value of < \cos6ij\ > obtained with the real 
data. Thus, the coincidence of these independent pre- 
ferred axes in such a small angular region is a highly 
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Fig. 6: A histogram of the distribution of < \cos6ij\ > as ob- 
tained from the Monte Carlo data, superpos ed w ith the value 
obtained from the real data of Table 2 (eq. (3.2)). Less than 
1% of the Monte Carlo data exceed the value corresponding 
to the real data. 



unlikely event. Even if we ignore the axes related to the 
CMB the coincidence of the three remaining axes in such 
a small angular region is still a relatively unlikely event 
with probability about 7% (Fig. 6). In fact for the three 
remaining axes of Table 2 we find 



< IcosOijl >=0.76 
while the Monte Carlo simulation gives 



< 



\cosOij\ >- 



0.5 ±0.16 



(3.4) 



(3.5) 



i.e. the real data are about 1.5a away from the Monte 
Carlo mean value. 

We therefore conclude that even though each of the 
axes of Table 2 does not by itself constitute statistically 
significant evidence for a cosmological anisotropy, their 
coexistence in a relatively small angular region is a very 
unlikely event which is most probably attributed to ei- 
ther an undiscovered physical effect or to a common basic 
systematic error that has so far escaped attention. 



4. CONCLUSIONS 

The main conclusions of our study may be summarized 
as follows: 

• The hemisphere of maximum accelerating expan- 
sion of the universe according to the Union2 data 



has a pole in the direction (/, b) = (309°lf , 18°lio) 
while the hemisphere of minimum acceleration is in 
the opposite direction. 
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Fig. 7: A histogram of the distribution of < \cosOij\ > as 
obtained from the Monte Carlo data, superposed with the 
value obtained from the real data of Table 2 excluding the 
axes related to the CMB (eq. (lO}). About 7% of the Monte 



Carlo data exceed the value corresponding to the real data. 



• The level of anisotropy of the Union2 dataset is 
larger than about 70% of simulated datasets but it 
is still consistent with statistical isotropy. 

• The coincidence of the anisotropy axes of Table 2 
in a relatively small angular region is statistically a 
highly unlikely event that hints towards a physical 
or systematic connection among the axes of Table 
2. 

The combination of these cosmological observations may 
involve some a posteriori reasoning in the sense that there 
are other observations that show no hint of a preferred 
axis. However, most large scale cosmological observa- 
tions may be analyzed in the context of a monopole, a 
dipole term (and higher moments) even if the significance 
level of a preferred direction is low and consistent with 
statistical isotropy. A real a posteriori reasoning would 
select some of these observations and points out that they 
point in similar directions. On the other hand our analy- 
sis includes all large scale cosmological observations that 
we could find in the literature which involve a preferred 
cosmological direction (even if it is of low significance by 
itself) and point out that these directions are abnormally 
close to each other. 

The confirmation of the existence of a cosmological 
preferred axis would constitute a breakthrough in cosmo- 
logical research. Given the present status of cosmological 
observations such a confirmation is one of the most prob- 
able directions from which new physics may emerge. 

Given the preliminary evidence for anisotropy dis- 
cussed in the present study, it is important to extend 
and intensify efforts for the possible confirmation of this 
evidence. Such confirmation may be achieved by extend- 
ing the Snia compilations towards larger datasets and 
deeper redshifts that span as uniformly as possible all 
directions in the sky. This is important in view of the 
fact that the Union2 compilation is less uniform and de- 
tailed in the south galactic hemisphere. In addition it 
is important to extend other cosmological data related 
to CMB low multipole moments, bulk velocity fiows and 
quasar polarization to confirm the present existing evi- 
dence for preferred axes in these datasets. Finally, alter- 
native probes of cosmological anisotropics may be consid- 
ered like higher CMB multipole moments, non-gaussian 
features and polarization in the CMB maps, alignments 
of geometric features of various structures on large scales 
(there is already some preliminary evidence for alignment 
of handedness of spiral galaxies [78 along an axis not far 
from the directions of the other preferred axes of Table 
2), alignment of optical polarization from various cosmo- 
logical sources or studies based on cosmic parallax [79 . It 
is also important to derive observational signatures that 
can clearly distinguish between the various different ori- 
gins of the preferred axes discussed in the introduction. 

Numerical Analysis Files: The Mathematica 7 files 
with the data used for the production of the figures along 
with a PowerPoint file with additional figures may be 
found at http://leandros.physics.uoi.gr/anisotropyi . 
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